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Abstract 



We formulate the data analysis problem for the detection of the New- 
tonian waveform from an inspiraling compact-binary by a network of arbi- 
trarily oriented and arbitrarily distributed laser interferometric gravitational 
wave detectors. We obtain for the first time the relation between the opti- 
mal statistic and the magnitude of the network correlation vector, which is 
constructed from the matched network-filter. This generalizes the calculation 
reported in an earlier work ( gr-qc/9906064 ), where the detectors are taken to 
be coincident. 
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I. INTRODUCTION 



Inspiraling compact binaries form prime candidates for detection by earth-based interfer- 
ometric gravitational- wave (GW) detectors owing to the well understood waveform (chirp) 
emitted by them. Searching for chirps using a network of such detectors is gaining impor- 
tance due to (a) its superior sensitivity vis a vis that of a constituent detector [1] and (b) 
improving feasibility for a real-time computational search. Here, we formulate the problem 
of how to optimally detect the Newtonian chirp using a network of arbitrarily orientated 
and arbitrarily located detectors. This extends a similar study in Ref. [2] of a network with 
coincident detectors. 

We use the maximum likelihood method for optimizing the detection problem. [3] A 
single likelihood ratio (LR) is deduced for the entire network. A super-threshold value 
for the maximized likelihood ratio (MLR) implies a detection. The MLR is obtained by 
maximizing the LR over the eight parameters that determine the Newtonian chirp: the 
distance r to the binary the inital phase 5 of the waveform, the polarization angle ip, the 
inclination angle e of the binary orbit, the time of arrival t a at a fiducial detector (fide), the 
source-direction angles {<p,6}, and the chirp time £. In principle, this can always be done 
numerically using a grid in the eight dimensional parameter space. In practice, however, 
such a strategy is computationally unfeasible and wasteful. We show that maximization of 
the LR over four parameters, {r, 5, ip, e}, can be performed analytically using the symmetries 
in detector responses. This allows us to scan this parameter subspace continuously. Further, 
the Fast Fourier Transform (FFT) can be used to maximize LR over t a , as in the case of a 
single detector. Such an analytic maximization and the FFT allow us to save substantially 
on computational costs. Numerical maximization is required over the remaining parameters, 
{(f), 0, £}, which we discuss in a future work. Here, we follow the convention laid out in Ref. 
[2]- 



II. THE SIGNAL 

There are four distinctly different reference frames of interest, associated with the source, 
wave, fide, and a representative detector in the network. Physical quantities in these frames 
are related by orthogonal transformations, O^, which are defined in terms of three sets 
of Euler angles that specify the orientation of one frame with respect to another. [4] Let 
x be an arbitrary three-dimensional real vector. Then, x wavc = Oi(tp, e, 0)x source , xgd e = 
2 (0, 0)x wavc , and x detector = 0^ x (a, (3, 7)x fide , Here, the source axes have been chosen in 
accordance with Ref. [5]. 

The wave tensor Wij associated with any source can be expanded in terms of the STF-2 
tensors 34n m an arbitrary frame as [2]: 

v>*(t) = M \(h + (t) - ih x (t)) r 2 n yz + (h + (t) + ih x (t)) T_ 2 n yz] , (2.i) 

V 15 L J 

where i, j denote spatial indices, and h + and h x are the two GW polarizations in the 
transverse-traceless gauge, as measured in some given frame. The expansion coefficients 
T ±2 n are the Gel'fand functions, [2] which depend on the Euler angles through which one 
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must rotate that frame into the frame in which w ij is being analyzed. The above form 
suggests the definitions, e % i = ^8%/ 15 T 2 n 3^2n and = ^8^/15 T^y^ for the left- and 
right-circular polarization tensors, respectively. They obey, e[ * — a%3 a%3 



e Ri e L,R e L,Rij — 1; an d 



e L, R e R, L ij 



0, in any frame. Thus 

w ij (t) = Re 



(h + (t)+ih„(t))e% 



= 2k Re 



R(t)4 



(2.2) 



where k = v^/ r ( U P to a normalization factor) and R = (h+(t) + ih x (t))/(2/c). For a chirp, 
we define R in the source frame. [5] Then h +jX are GW amplitudes for a face-on binary (i.e., 
for e = 0), and R depends only on {5, t a ,£}. 

The response amplitude (i.e., the signal) in the 7-th detector is the scalar product s 1 = 
^djj, which depends on projections of e l [ R onto the J-th detector tensor, c? 7 -. One such 



w 



projection defines the extended beam-pattern function: 

F J = ei < $ j =Tf{i/,,e,0)D I p , p = ±2 
which corresponds to the left-circular polarization. Above, 
/8tt 



Di = V 15 V(0> 0)4^„ = ig 1 V(0, 61, 0) (l£ - T!_* 2n 



(2.3) 



(2.4) 



where T^ 2n 



T±2n(a/,/5/,7/) and gZ 7 - = ^ I (n{ i n{ j - n^n^), 



; . with nf 2 being unit vectors 

along the two arms of the J-th interferometer, respectively. Also, g 1 is the detector's noise 
power spectral density. [6] Then, 



s 7 (t) = 2k Re (F 7 *^ 7 ) = 2k Re (F J *S J e w ) 
where R 1 is defined via Eq. (|2.2| ) and S 1 is independent of 5. 



(2.5) 



III. THE OPTIMAL NETWORK STATISTIC 

Under the Ney man- Pearson decision criterion, [3] the optimal network statistic is the 
network LR, A. If the noise in each detector is additive and independent of the noise in any 
other detector in the network, then A reduces to a product of the individual detector LR's. [2] 
Further, for Gaussian noise, [7] the logarithmic likelihood ratio (LLR), In A, simplifies to the 
following sum of LLR's of N individual detectors [2]: 

N ^ N N ^ 

\ n \ = J2(si,X I ) I --J2( s Ii s l)l = b J2( Z Ii X l)l~n b2 » (3- 1 ) 
1=1 1 1=1 1=1 1 

where b = 2k(^Ji =1 ||-^i"|| 2 ) 1//2 and zi = s//b. Above, r appears only in b. 

2 

Maximizing In A with respect to b and S gives, lnA^ = J2f=iQiCj /2, where 
Qi = 2/tFj/b and C} = (Sj,xj)i.[2] This shows that the network vector S, with SV's as 
its components, is the matched network-filter. Also, lnA^ is a function of six parameters, 
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namely, {if}, e, t a , (f>, 9, £}. To extend these results to the case where the detectors are arbi- 
trarily located, note that the dependence of In A| ^ g on {if>, e} can be isolated. This is because 
the network vector Q, with Q 7 's as its components, is: 

Q = II F (T 2 - 2 (ifj, e, 0)D_ 2 + T 2 (if}, e, 0)D 2 ) = Q- 2 D_ 2 + Q 2 D 2 , (3.2) 

where D p define network vectors with D p as their components; D p are their normalized 

counterparts. Thus, Q2 = D2 • Q — Q +2 + Q 2 D 2 • D_ 2 . Hence, {D 2 , D_ 2 } define a two- 
dimensional complex plane, V (the helicity space, a subspace of C N ), on which a metric g pg 
can be defined. Then, Q p = g pq Q q with p,q — ±2. The iV-dimensional correlation vector C, 
in general, lies outside V. However, Q lies totally in V. Thus, the statistic reduces to 

lnA| M = |Q-C| 2 /2= |Q-q,| 2 /2 , (3.3) 

where C-p is the projection of C onto the helicity space. 

Maximization of In s over {ip, e} is achieved by aligning Q along C-p. This requires that 
Q ±2 = Cj5 2 /\\ C-p \\ v , which implies: Q +2 /Q~ 2 = C^ 2 jC v 2 . Since the RHS above can take 
any value in the complex plane, we need to prove that Q +2 /Q~ 2 obeys the same property. 
To this end, note that, 



Q+ 2 T 2 " 2 (^,e,0) /1-cose^ 2 



:) exp(4#) , (3.4) 



Q- 2 T 2 2 (^,e,0) VI + cose, 

which can indeed attain any value on the Argand plane. The values of if) and e that maximize 
the statistic are, ip = arg(x)/4 and e = cos _1 [(l — \J\\x\\ )/(l + \J\\x\\ )] where x = C^ 2 /C^ 2 . 
Thus, the statistic maximized over these four parameters is, 

lnA| M ^=||C P || 2 /2. (3.5) 

Let V* 1 denote a pair of orthonormal, complex basis vectors on V. Then, 

II Cp || 2 = || C+ || 2 + I) C~ || 2 = (c +) 2 + (c+ /2 ) 2 + (c -) 2 + (c~ /2 ) 2 , (3.6) 

where C ± = C-p ■ V ± = c^ + ic^, 2 . It can be verified that the network statistic is, therefore, a 
sum of the squares of four Gaussian random variables with constant variance. This simplifies 
the computation of thresholds and detection probabilities. 

A network filter is constructed as follows: For a given £ compute the Newtonian chirp 
for t a = 0. Then, for a given direction {0, 9}, use the appropriate time-delays with respect 
to fide to time-displace the chirp at each detector. This collection of time-displaced chirps 
constitute the network filter. Also, t a is obtained by shifting the network filter 'rigidly' on 
the time axis, which can be done efficiently using FFT. The bank of filters on {0, 9,1;} can 
be obtained by correlating two neighboring normalized filters in the usual way. This work 
is now in progress. 
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IV. CONCLUSION 

We have given here a formulation to optimally detect the Newtonian chirp with a network 
of detectors in which the noise is additive, Gaussian, and uncorrelated between detectors. 
We have shown how this can be done efficiently by analytically maximizing the LR over 
four parameters and using FFT to maximize over the time-of-arrival. In a future work we 
hope to address key issues such as required computational power for such a search and also 
estimate errors in parameter values. 
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